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ABSTRACT 

ln this paper, damage ¡n simple and double-lined tunnels is modelled as the 
occurrence of hinges. Primary and secondary linings are modelled with a beam- 
column element with embedded discontinuities, which was developed by Juárez- 
Luna y Tenorio-Montero (2019). In this element, damage ¡s modelled as the 
development of hinges under flexural forces, transverse dislocations under shear 
torces and axial discontinuities under axial forces. This element uses constitutive 
models of reinforced concrete prismatic sections to inelude independently the 
capacities of the moment, shear forcé and axial forcé, based on damage mechanics. 
This beam-column element with embedded discontinuities was calibrated with an 
experimental test of a simple tunnel lined reported in the literature, the numerically 
computed load-displacement curve was in agreement with that of the experimental 
test. Numerical examples of simple and double-lined tunnels were carried out, 
showing the capability of this finite element for modelling damage. In these numerical 
examples, ovalisation load vs ovalisation curves were computed, showing that peak 
ovalisation load-carrying capacity in the double-lined tunnel ¡s approximately the 
double than the peak ovalisation load-carrying capacity in the single-lined tunnel, but 
the latter is more flexible. 
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MODELADO DEL DAÑO EN TÚNELES DOVELADOS 
SIMPLES Y CON DOBLE REVESTIMIENTO MEDIANTE 
ARTICULACIONES EMBEBIDAS EN ELEMENTOS 
GRUESOS VIGA-COLUMNA 

RESUMEN 

En este artículo se modela el daño en túneles con revestimiento simple y 
doble como la ocurrencia de articulaciones. El revestimiento primario y secundario 
se modela con el elemento viga-columna con discontinuidades interiores, 
desarrollado por Juárez-Luna y Tenorio-Montero (2019). En este elemento, el daño 
se modela como el desarrollo de articulaciones bajo fuerzas flexionantes, 
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dislocaciones bajo fuerzas cortantes y discontinuidades axiales bajo fuerzas axiales. 
Este elemento utiliza modelos constitutivos de secciones prismáticas de concreto 
reforzado para incluir independientemente las capacidades de momento, fuerza 
cortante y fuerza axial, con base en la mecánica del daño. Este elemento viga- 
columna con discontinuidades interiores se calibró con una prueba experimental de 
un túnel con revestimiento simple reportada en la literatura, en el que la curva 
calculada de la carga contra el desplazamiento es congruente con aquella de la 
prueba experimental. Ejemplos numéricos de túneles con revestimiento simple y 
doble se realizaron, mostrado la capacidad de este elemento finito para modelar el 
daño. En estos ejemplos se calcularon las curvas de cargas de ovalización contra 
la ovalización, mostrando que la capacidad de carga máxima en el túnel con doble 
revestimiento es aproximadamente el doble que la capacidad de carga máxima que 
en el túnel con revestimiento simple, pero el último es más flexible. 


Palabras clave: Túnel; daño; articulación; dislocación; discontinuidades 

1 Introduction 

El daño en túneles se atribuye a cargas actuantes, geometría de la estructura, 
comportamiento constitutivo de los materiales, entre otros. El incremento del daño 
en un túnel puede producir su colapso. El daño en túneles ocurre en distintas 
condiciones, e.g., durante cargas sísmicas, estos se someten a cargas axiales y 
cortantes que resultan en distintos mecanismos. Estos mecanismos de colapso 
dependen de la geometría y de su revestimiento, pues túneles continuos y 
segmentados exhiben mecanismos de falla distintos. El estudio de la ocurrencia del 
daño y su evolución en túneles es importante para conocer la capacidad máxima y 
la capacidad residual. 

El método de los elementos finitos (MEF) es probablemente uno de los métodos 
computacionales más utilizados para modelar la ocurrencia del daño y su evolución 
en elementos estructurales, donde éstos se modelan detalladamente con elementos 
sólidos o con modelos más simples mediante elementos viga-columna. Estos 
elementos estructurales se modelan con elementos finitos (EF) a los que se les 
asignan modelos constitutivos, los cuales incluyen comportamiento no-lineal. 
Cuando una estructura se modela con EF sólidos 3D, el daño se modela como la 
ocurrencia de grietas o la degradación del material, pero si la estructura se modela 
con EF viga-columna ID, el daño se modela como como la ocurrencia de 
articulaciones. Con esta última aproximación se tienen resultados razonables y se 
reduce el tiempo computaclonal, lo cual se atribuye a que se requieren menos 
grados de libertad para modelar una estructura que modelarla con EF sólidos; 
además, el comportamiento constitutivo de las secciones transversales de concreto 
reforzado se simplifica al modelo constitutivo de una articulación. 

Para modelar el daño en estructura, Baker y Hyman (1969) modelaron la ocurrencia 
del daño como articulaciones en elementos viga-columna, utilizando relaciones 
constitutivas momento-curvatura con un comportamiento elasto-plástlco con 
endurecimiento Isótropo. Sin embargo, esta formulación no incluye el ablandamiento 
después de alcanzar la carga máxima. Ehrlich y Armero (2004) formularon un EF 
viga de Timoshenko con discontinuidades embebidas, el cual es capaz de modelar 
saltos de rotación y transversales; sin embargo, este elemento necesita operadores 
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para estabilizar su solución. Este EF proporciona una matriz de rigideces no 
simétrica, la cual puede tener estabilidades numéricas. Estos mismos autores, 
Armero y Ehrllch (2006), desarrollaron un EF con discontinuidades embebidas 
basado en la teoría de vigas de Euler-Bernoulli; este elemento tiene la capacidad 
de modelar sólo saltos de rotación. 

Dujc et al. (2010) desarrollaron un modelo multlescala para estudiar 
discontinuidades embebidas axiales y rotacionales, utilizando EF en 2D para el 
modelo en mlcroescala y EF vigas de Euler-Bernoulli para el modelo macroescala. 
Juárez y Ayala (2012) desarrollaron elementos vigas gruesos y delgados con 
discontinuidades embebidas con base en funcionales de energía, cuya 
aproximación con el MEF proporciona matrices de rigideces simétricas. Estos 
elementos finitos con discontinuidades interiores fueron validados con resultados 
experimentales y numéricos reportados en la literatura. Estos autores además 
desarrollaron modelos constitutivos momento flexionante contra salto de rotación y 
fuerza cortante contra salto de desplazamiento transversal, el cual incluye 
ablandamiento después de alcanzar su valor umbral. Juklc et al. (2013) 
desarrollaron un elemento finito con discontinuidades embebidas con base en la 
teoría de vigas de Euler-Bernoulli, el cual modela sólo el salto de rotación. La 
relación constitutiva momento flexionante contra curvatura es similar al modelo 
desarrollado por Dujc et al. (2010), el cual modela el ablandamiento después de 
alcanzar el momento último. Juklc (2013) desarrolló un EF, Incluyendo el 
comportamiento no lineal axial en vigas de Tlmoshenko y de Euler-Bernoulli. La 
sección transversal de la viga se modeló con fibras, las cuales tienen un modelo 
constitutivo esfuerzo contra deformación. Las fibras presentan esfuerzos axiales 
producidos por los momentos flexlonantes. 

Los túneles son un tipo particular de estructuras con forma ovalada. Cuando los 
túneles se construyen con máquina tuneladora, éstos se construyen principalmente 
con dovelas prefabricadas como revestimiento primarlo, y algunos de ellos tienen 
además un revestimiento secundarlo de concreto colado en sitio. Blom (2002) trata 
con el diseño del revestimiento dovelado de túneles en suelos suaves. Él desarrolló 
una aproximación analítica para calcular los desplazamientos radiales y el momento 
flexionante tangencial en los anillos, el cual incluye la rigidez rotacional de las juntas 
longitudinales y la Interacción lateral entre anillos para un revestimiento en un medio 
continuo elástico. Los resultados del desplazamiento radial y del esfuerzo radial 
tangencia calculados con la solución analítica, modelado con EF viga-columna y con 
EF sólidos son congruentes con los valores medidos en la prueba a escala natural 
reportada por Blom y van Oosterhout (2001). Wlnkler etal. (2004) realizaron análisis 
no-lineales de modelos de EF del revestimiento de un túnel construido con dovelas 
hexagonales prefabricadas. Estos autores mejoraron un modelo de plasticidad para 
concreto simple, acoplando daño en tensión con daño en compresión e 
Introduciendo un modelo de daño escalar Isótropo. Las trayectorias de 
agrietamiento calculadas muestran congruencia con el agrietamiento observado en 
sitio. Luttlkholt (2007) además desarrolló modelos de EF de los especímenes a 
escala natural reportados por Blom y van Oosterhout (2001). Él modeló los 
especímenes de los túneles con EF en 1D, utilizando EF tipo viga y modelos en 2D, 
utilizando elementos de esfuerzo plano. Arnau y Molins (2011) realizaron 
simulaciones numéricas de una prueba en sitio. Esta prueba fue modelada con dos 
aproximaciones distintas: un modelo es esfuerzo plano en 2D y con elementos 
cascarón en el programa Diana. En el modelado de la respuesta estructural del 
revestimiento se Incluyeron cuidadosamente el comportamiento post-agrletamlento, 
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juntas longitudinales y la interacción suelo-estructura; además, se modelaron las 
fibras de acero como refuerzo en el concreto. Los desplazamientos calculados 
numéricamente en el revestimiento coincidiendo adecuadamente con las medidas 
experimentales, los cuales fueron producidos por las juntas y por la rotación de 
secciones agrietadas. Las similitudes entre las trayectorias de agrietamiento 
obtenidas numéricamente con las reportadas en sitio corroboran la conveniencia de 
los modelos utilizados, van der Waart van Gulik (2012) desarrollaron análisis no 
lineales de modelos de anillos segmentados circundados por suelo, concluyendo 
que los túneles dovelados con mayor diámetro son más sensibles a inestabilidad 
que los túneles con diámetros menores. Él modeló el anillo de un túnel con EF tipo 
viga y las juntas longitudinales con resortes rotacionales. Arnau y Molins (2015) 
estudiaron el mecanismo de interacción de la sección real de un túnel en escenarios 
múltiples. Ellos realizaron modelos numéricos con elementos finitos en 3D para 
simular la respuesta estructural no lineal de once anillos en el programa Diana 9.4. 
Las dovelas se modelaron con EF tipo cuadriláteros cascarones curvos con 4 nodos, 
las juntas longitudinales con elementos de interfaz para elementos cascaron, la 
interacción suelo-estructura con elementos colocados en las direcciones radial, 
tangencial y longitudinal. Se concluye que la rigidez del suelo circundante determina 
el grado de interacción posible entre los anillos adyacentes, mientras que la 
compresión longitudinal del revestimiento determina la carga máxima para la cual 
se produce una Interacción completa. Galván etal. (2017) investigaron la respuesta 
estructural del revestimiento de túneles dovelados al avance de una máquina 
tuneladora durante la construcción de un túnel típico en México. Ellos realizaron un 
estudio paramétrico considerando el comportamiento no-lineal de los materiales, la 
Influencia del número de anillos, la carga axial residual y la posición de las juntas 
continuas. Estos autores concluyen que la carga residual incrementa la capacidad 
estructural del revestimiento. Además, se muestra que el mecanismo de colapso se 
modifica con el acoplamiento entre los anillos. Los análisis numéricos de los 
modelos de túneles los realizaron en el programa Ansys. Las dovelas se modelaron 
con EF tipo sólido, pero las juntas longitudinales y transversales se modelaron con 
elementos de contacto. Recientemente, Juárez-Luna y Tenorio-Montero (2019) 
desarrollaron un elemento finito con discontinuidades embebidas, el cual modela 
discontinuidades de desplazamiento axial, rotación y desplazamiento transversal. 
Este elemento fue desarrollado con base en funcionales de energía de una barra y 
de una viga gruesa, cuya aproximación proporciona una matriz de rigideces 
simétrica. El acoplamiento de los grados de libertad de elementos tipo barra y tipo 
viga proporciona un elemento viga-columna que Incluye desplazamiento axial, 
rotación y desplazamiento transversal, así como sus respectivos saltos. Estos 
autores mostraron la capacidad de este elemento finito al modelado del daño de 
túneles dovelados con revestimiento simple y revestimiento doble, en la que el daño 
evoluciona como la ocurrencia de articulaciones. 

El contenido de este artículo es como se describe a continuación. En la sección 2, se 
presentan las matrices de rigideces tangentes del elemento viga-columna con 
discontinuidades interiores desarrollados por Juárez-Luna y Tenorio-montero (2019). 
En la sección 3 se proporcionan los modelos constitutivos para describir el 
comportamiento de las dovelas: las fuerzas contra los desplazamientos continuos y 
las fuerzas contra los saltos de los desplazamientos. Además se proporcionan el 
modelo constitutivo del momento contra la curvatura en las jutas longitudinales. En la 
sección 4 se muestran ejemplos numéricos de túneles con revestimiento simple y con 
doble revestimiento que muestran la aplicabilidad de este elemento para modelar el 
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daño. Finalmente, en la sección 6, se presentan las conclusiones derivadas de este 
trabajo. 

2 Elemento viga-columna con discontinuidades embebidas 

El elemento finito viga-columna con discontinuidades embebidas fue 
desarrollado por Juárez-Luna y Tenorio Montero (2019). Este elemento tiene en total 
nueve grados de libertad (GDL), tres por nodo; los GDL de ambos extremos están 
asociados a los desplazamientos axiales, u, a los desplazamientos transversales, v, 
y a los giros 0. Por otra parte los GDL del nodo Interno está asociado a los saltos 
de desplazamiento axial, [|u|], trasversal, [|tv|], y giro, [|0|], donde existe 
concentración de deformaciones, como se muestra en la Figura 1. 


sj 2 

, n 


Q~ S fi + b) - 

Figura 1. Elemento viga columna con discontinuidades embebidas: a) GDL y b) 

desplazamientos 


La formulación de Juárez-Luna y Tenorio-Montero (2019) la matriz de rigideces del 
elemento viga-columna gruesa con discontinuidades Interiores es: 



Donde: 
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-hW\) 

( 2 ) 


Donde: 

Cj=av([ M ])/a[ H ] y Cj = 5A/([|0|])/d[|0|] son los operadores tangentes de 
salto de desplazamiento transversal y de salto de rotación, respectivamente. 

En el elemento con discontinuidades Interiores, los GDL en el nodo Interno (ver 
Figura 1) pueden condensarse como grados de libertad secundarios, y expresarlos 
en términos de los GDL (d, w,6) en ambos extremos del elemento. La matriz de 
rigideces en la ecu. (1) se reescribe utilizando una partición de submatrlces tal que: 



Donde: 


K ee , K e [| e |], K¡ie|]e y K¡ie|][|e|] son submatrlces asociadas a sus correspondiente GDL, 
incluyendo las discontinuidades o saltos en el elemento. F n y ñ[| e |] son los vectores 
que contienen las fuerzas nodales y momentos. A d e son los desplazamientos 
nodales y A[|e|] son los saltos en la discontinuidad. Estas matrices particionadas 
son: 
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En la ec. (3), los GDL pueden condensarse, de acuerdo con Tena-Colunga (2007), 
como: 


F, = K„M r + K t [k|] A[>|] 

(5) 

%o = í: [i# íW ' +a: [h][h) a D <, I] 

(6) 


Resolviendo para A[|e|] en ec. (6), se tiene: 

A [H] = K mñ\i ~ K ím K me Ade (7) 

Substituyendo ec. (7) en ec. (5), se tiene: 

K ee ~ K e[\e\] K me\] K Me ] l Ad e ] = [ F n ~ VoÍm^HI ] (8) 

La cual se representa como: 

WM=[*1 o» 

La matriz de rigideces condensada, Ka, dada en la ec. (8), es: 

l*\l 

Y el vector condensado que contiene las fuerzas residuales R* es: 


(11) 
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La matriz de rigideces del elemento viga-columna con discontinuidades interiores, 
dado en la ec. (1), puede condensarse estáticamente con base en la ec. (8). 
Consecuentemente, solo los GDL en los extremos del elemento se activan, pero los 
GDL de los saltos se condensan de acuerdo a las ecs. (9) a (11). 

3 Modelos constitutivos 

Los modelos constitutivos representan el comportamiento idealizado de los 
materiales, incluyendo su intervalo elástico y no-lineal. En los modelos de daño 
continuo, la degradación del material se distribuye en el volumen del elemento finito 
dañado. Por el contrario, en un modelo de daño discreto, la degradación del material 
se concentra en una zona de espesor cero dentro de un elemento finito. 

En el caso de elementos finitos tipo viga, se requieren relaciones constitutivas 
continuas entre las fuerzas y los desplazamientos o relaciones constitutivas 
discretas entre las fuerzas y los saltos de desplazamientos. Un ejemplo de un 
modelo de daño continuo ¡sotrópico para elementos tipos vigas son los formulados 
por Juárez y Ayala (2012), dados por las ecuaciones siguientes: 


Densidad de energía libre 
Ecuación constitutiva 


y/(X,r) = (\-d(r))y 0 

8x 


Variable de daño 


Ley de evolución 


Criterio de daño 

Ley de endurecimiento 
Condiciones de carga 
y descarga 


d(r) = 1 - —; tfe[0,lj; ge[0,r o ] 
re[r 0 ,co] 

i 0 «=» Ve 

f{r F ,q) = T F -q = ^F:C~ l :F -q\ j‘ 7 ‘ E [°’ r °] 

[ C i | f-0 r o 

4^H d (r)& H d (r) = q\r) < 0 
f(T F ,q)< 0; r>0; yf(r F ,q) = 0 

rfir F ’<l) = 0 

( 12 ) 


donde V es la densidad de energía libre, el área bajo la curva del modelo 
constitutivo, / es el vector de desplazamiento, Cson las constantes elásticas, Fes 
el vector de fuerza, de s la variable de daño, la cual es dependiente de la variable 
de ablandamiento q, He s el parámetro de ablandamiento y y es el multiplicador de 

daño que determinan las condiciones de carga-descarga. La función f(zF,q) delimita 
el dominio elástico definiendo la superficie de daño en el espacio de las fuerzas. 
Por otra parte, el modelo de daño discreto para vigas es: 
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Densidad de energía libre ^([|z<|],a) = {\-co)[// 0 ([| M |])> j^" 0 ítl M Ü) lo ^'([l“l]) í ^[M] 

dy/([|w|],a) 


Ecuación constitutiva 

Variable de daño 

Ley de evolución 

Criterio de daño 
Ley de endurecimiento 
Condiciones de 
carga-descarga 


F = - 


5[|u|] 


- = (l — ®)C |m|] 


i ?(«) r ii 
tu = l-——; oe[-oo,lj 


t$=/l=— (a); ae(0,oo) 

f(T,q) = r T -q\ v T = ||f||(C) ‘ = <Jf[C]- 1 F 
${a) = H$; H = q\a) < O 

/ < 0; X > 0; X/ = 0 
X/=0 

(13) 


donde co es la variable de daño, la cual depende de la variable de ablandamiento 
discreto q. El resto de las variables se definen con el mismo nombre que el modelo 

de daño continuo, pero Incluyendo el término discreto, e.g., H es el parámetro de 
ablandamiento discreto. 


La razón de cambio de las fuerzas con respecto a los desplazamientos corresponde 
al operador constitutivo tangente C T , de un modelo constitutivo continuo. Por otra 
parte, la razón de cambio de las fuerzas con respecto a los saltos de 
desplazamientos corresponde al operador constitutivo discreto, CJ, de un modelo 
constitutivo discreto. Estos operadores constitutivos tangentes respectivamente 
son: 


r T = dF r T = dF 

dz’ 


(14) 


Cuando el material está en condición de carga en el Intervalo no-lineal, el operador 
constitutivo tangente se calcula con la ec. (15), pero para el intervalo elástico o en 
condición de descarga, estos operadores tangentes de calculan con la ec. (16). 


C T =(1 -d)C- C ±-p-{C: X ®X-C), C' d =(1 -a))C-^^(c{\u\]®[\u\}c) 

(15) 

C T =(l-d)C, Cj =(l-<s>)C (16) 

El elemento viga-columna requiere de un modelo constitutivo que relaciona la fuerza 
normal, N, la fuerza cortante, V, y el momento flexlonante, M, con la deformación 
normal, e, deformación por cortante, y, y la curvatura, k, respectivamente. Para el 
caso de los elementos viga-columna con discontinuidades desarrolladas por Juárez- 
Luna y Tenorio Montero (2019), estos requieren de modelos de daño discreto, N- 
[|t/|], l/-[|i/v|], M-[|0|], Estos modelos constitutivos en sus versiones continua y discreta 
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se muestran respectivamente en las Figuras 2, 3 y 4, donde los subíndices yy use 
refieren a las fuerzas de fluencia y última, respectivamente. 


Níe) N ([IkI]) 



Figura 2. a) Relaciones constitutivas continua fuerza normal contra deformación 
normal y b) relación constitutiva discreta de la fuerza normal contra el salto de 

desplazamiento axial. 


V([lwl]) 




Figura 3. a) Relación constitutiva continua fuerza cortante contra deformación 
cortante y b) relación constitutiva discreta de la fuerza cortante contra el salto de 

desplazamiento transversal. 


M( K) 



a) 


MdieiD 



b) 


Figura 4. a) Relación constitutiva continúa del momento flexlonante contra la 
curvatura y b) relación constitutiva discreta del momento flexlonante contra el salto 

de rotación. 
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El comportamiento constitutivo de la junta longitudinal en túneles se describe con 
los resultados del momento contra la rotación reportados por Leonhardt y Relmann 
(1966), quienes probaron dos piezas de concreto con una junta como se muestra 
en la Figura 5a. La junta tiene un ancho a y una longitud b. La formulación 
matemática de este modelo constitutivo fue desarrollado por Jansen (1983), que 
describe la curva del momento contra la rotación, mostrada en la Figura 5b, con los 
valores dados en la ec. (17). 



Figura 5. a) Junta probada por Leonhardt y Relmann (1966) y b) curva momento 
contra curvatura (Janssen, 1983). 


C = 


bcrE„ 


9ba\ 


hC = ■ 


12 
2M 
^ Na 


-1 


8N 


0 < 


2 N 
Eba 


Ai; 0 > 


2N 

E.ba 


(17) 


donde E c es el modulo elástico del concreto, M es el momento flexlonante y N es la 
fuerza normal. El comportamiento constitutivo de la junta longitudinal se Idealiza con 
la curva bilineal mostrada en la Figura 6, en la que el momento se calcula como: 


M ={\-d)C r 6 


(18) 


M(0) 



Figura 6. Curva bilineal del momento vs la rotación. 
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El modelo constitutivo con discontinuidades interiores y los modelos constitutivos se 
implantaron en subrutinas del programa FEAP (Finite Element Analysis Program) 
desarrollado por Taylor (2008). Estas subrutinas se escribieron en el lenguaje de 
programación Fortran. 


4 Casos de estudio 


4.1 Túnel con revestimiento simple: Túnel Botlek Railway 

Este caso de estudio consiste de un espécimen a escala natural de tres 
anillos del túnel Botlek Railway, el cual fue ensayado en el Laboratorio Stevln II de 
la Universidad de Tecnología de Delft por Blom y van Oosterhout (2001). El 
espécimen tiene tres anillos con un ancho de 1.5 m formados por dovelas con 
espesor de 40 cm, como se muestra en la Figura 7a. Cada anillo tiene un diámetro 
de 8.65m, que está compuesto por siete segmentos y una dovela llave situada en 
posiciones distintas en cada anillo. El anillo medio, identificado como anillo 1, tiene 
su primera junta longitudinal en la posición 0.5 y la dovela de remate en la posición 
6.5. Los anillos superior e inferior, identificados como anillo 2, tiene la primera junta 
longitudinal en la posición 0 y la dovela llave en la posición 5.5 (Blom 2002), como 
se muestra en Figura 7b. Las juntas entre las dovelas están alineadas entre el anillo 
superior e inferior, anillos 2. En las juntas entre anillos se utilizó madera 
contrachapada como material de empaque. Por otra parte, en las juntas 
longitudinales no se utilizó material de empaque, teniendo una superficie de 
contacto tipo concreto-concreto. 



Figura 7. a) espécimen a escala natural y b) 
en geometría (adaptada 


Localización de la primera 
junta longitudinal del anillo 2 



localización de juntas longitudinales 
de Blom 2002). 


El espécimen a escala real se cargó para dos distintos valores de fuerza axial. Una 
carga axial alta con interacción entre los tres anillos, experimento COI, y una carga 
axial con baja interacción entre los tres anillos, experimento C02. Considerando la 
baja interacción de este último experimento, se realizó el modelado del anillo central, 
anillo 2, como una aproximación a la baja interacción de los anillos. De acuerdo a la 
prueba C02, se aplicó primeramente una carga radial uniforme de 654.94 kN por 
actuador hidráulico. Posteriormente, se aplicó una carga de ovalización de 25 kN 
por actuador, resultando las fuerzas por actuador mostradas en la Figura 8. La carga 
uniforme produce fuerzas normales, contracción del anillo, pero manteniendo la 
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forma circular. Por otra parte, las cargas de ovallzaclón producen momentos 
flexlonantes, que flexlonan al anillo y que produce una forma ovalada del anillo. 




Figura 8. Anillo con: a) actuadores radiales y b) cargas máximas aplicadas (kN). 

Cada dovela tiene como refuerzo longitudinal trece barras de refuerzo con un 
diámetro de lOmm (10.21 cm 2 ) en la parte superior e Inferior, como se muestra en 
la Figura 9a. El refuerzo transversal consiste de estribos con un área total de 2.85 
cm 2 . La junta longitudinal tiene un ancho de b=1.39 m y un largo a=0.17 m, como 
se muestra en la Figura 9b. Las propiedades mecánicas del concreto son: módulo 
de Young de E c =40 GPa y esfuerzo último a compresión f’ c =27 MPa. El acero de 
refuerzo tiene las propiedades mecánicas siguientes: módulo de Young E s =200 GPa 
y esfuerzo de fluencia de f y =435 MPa. 


b) 



Figura 9. Propiedades: a) sección transversal de la dovela y b) junta longitudinal. 


La malla del anillo del túnel se muestra en la Figura 10, cada dovela se mallo con 
ocho elementos finitos tipo viga-columna (56 elementos con las propiedades 
mecánicas de las dovelas) en el programa GID. Las juntas entre las dovelas se 
modelaron con elementos tipo resorte rotacional (elementos número 57 a 63), a los 
cuales se les asignó el comportamiento constitutivo de Jansen (1983). La dovela 
llave se modeló con un solo elemento ente las dovelas vecinas, al cual se le asignó 
el elemento número 58. Las relaciones constitutivas asignadas a las dovelas y a las 
juntas longitudinales se proporcionan en el apéndice Al. 
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Figura 10. Modelo numérico en 2D del anillo 2 del túnel Botlek Rallway. 

La curva de la carga de ovallzaclón contra el desplazamiento calculada con los EF 
con discontinuidades embebidas muestran concordancia con la curva experimental 
reportada por Blom y van Oosterhout (2001), como se muestra en la Figura 11 a. La 
primera articulación ocurre en la parte inferior a una carga de 13.78 kN/actuador (1). 
Posteriormente, dos articulaciones ocurren en la parte superior a una carga de 18.20 
kN/actuador (2) y finalmente, dos articulaciones ocurren en ambos lados del túnel a 
una carga de 23.43 kN/actuador (3) como se muestra en la Figura 11b. Se observa 
que la ocurrencia del daño modelado como articulaciones con elementos viga- 
columna con discontinuidades embebidas son congruentes con las grietas 
calculadas con el modelo de agrietamiento distribuido reportadas por Luttikholt 
(2007), pues las articulaciones y grietas ocurren en la misma posición dentro de las 
dovelas y de las juntas longitudinales, como se muestran en la Figura 11c. Los 
elementos mecánicos se muestran en la 

Figura 12, donde se muestra que las juntas longitudinales no alcanzan los valores 
umbrales de fuerza cortante y fuerza normal, 1920 kN y 6814.33 kN, 
respectivamente. Sin embargo, los momentos flexionantes calculados en las 
dovelas alcanzan los valores umbrales de 184.89 kN-m. 



a ) 

Figura 11. Curva de la carga de ovallzaclón contra el desplazamiento vertical, b) 
desarrollo de articulaciones en el túnel con revestimiento simple y c) comparación 
de articulaciones con patrones de agrietamiento ocurridos durante la prueba 
experimenta en el anillo 2 (adaptado de Luttikholt 2007). 
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a) b) 




Figura 12. Elementos mecánicos: a) momento flexlonante (kN-m), b) fuerza 
cortante (kN) y c) fuerza normal (kN). 

4.2 Túnel con revestimiento simple y doble 

Este modelo estudia el espécimen a escala real de un túnel dovelado con 
doble revestimiento, el cual fue probado en México por Agullar et al. (2011), el 
revestimiento primarlo fue construido de develas prefabricadas de concreto y el 
revestimiento secundarlo fue construido con concreto colado en sitio. El 
revestimiento primarlo de dovelas tiene un anillo medio, Identificado como anillo 1, 
con un ancho completo y otros dos anillos, Identificado como anillo 2, con la mitad 
del ancho en la parte superior e Inferior. El revestimiento secundarlo consiste de un 
anillo Interior colado en sitio que tiene un espesor constante en el centro y en la 
parte superior del túnel, pero tiene un espesor variable en la parte Inferior como se 
muestra en la Figura 13a. Este túnel tiene siete dovelas y una dovela clave como se 
muestra en la Figura 13b. Las dovelas tienen una longitud de 4.0 m y la dovela clave 
una de 1.86 m; todas las dovelas tienen un alto de 1.50 m y un ancho de 0.40 m. 

Cada dovela tiene un área de acero de refuerzo longitudinal de Ai=12.5 cm 2 en la 
parte superior e Inferior; además, éstas tienen un área de acero de refuerzo 
transversal de At=2.85 cm 2 como se muestra en la Figura 13c. El diámetro interno 
de los anillos es de d=9.11m, las juntas longitudinales tiene un espesor b=1.39 m y 
una altura a=0.17 m como se muestra en la Figura 9b. Las propiedades mecánicas 
del concreto son: módulo de Young E c =21 GPa, radio de Polsson u c =0.2 y esfuerzo 
último a compresión o m =23.5 MPa. Las propiedades mecánicas del acero de 
refuerzo son: módulo de Young E s =200 GPa, radio de Polsson u s =0.3 y esfuerzo de 
fluencia cr y =412 MPa. 
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Figura 13. Túnel con doble revestimiento: a) esquema, b) geometría y c) sección 
transversal de dovela (adaptado de Aguilar et al. 2011). 


Este espécimen fue probado en dos etapas. Durante la primera etapa, el espécimen 
sólo tenía el revestimiento primarlo de dovelas sujeto a las cargas mostradas en la 
Figura 14a. En la segunda etapa, el revestimiento secundario de concreto colado 
en sitio, el cual tiene un espesor de 0.20 m a los lados y en la parte superior, pero 
tienen un espesor de 1.1 m en la parte más ancha en la parte inferior del anillo, 
como se muestra en la Figura 14b. Las cargas aplicadas tienen un coeficiente 
lateral del suelo de 0.83 en la primera etapa y de 0.61 en la segunda etapa; estas 
cargas se aplicaron mediante doce actuadores radiales. Las dovelas fueron 
sujetadas a un muro de reacción circular interior mediante unas barras pretensadas 
radlalmente, las cuales se colocaron a cada 30°. 

a) b) 

740.46 1005.5 




Figura 14. Esquema de cargas durante prueba de túnel: a) en la primera etapa 
(sólo revestimiento externo) y b) en la segunda etapa (revestimiento interno y 
externo) (adaptado de Aguilar et al. 2011) 

En la prueba experimental, doce puntos de control se asignaron como se muestra 
en la Figura 15a, teniendo la configuración deformada mostrada en la Figura 15b, 
donde se tiene una acortamiento entre los puntos uno y siete en la dirección vertical, 
pero se tiene un alargamiento entre los puntos 10 y 4 en dirección horizontal. Los 
resultados numéricos de las curvas de la carga vertical contra el acortamiento y la 
curva de la carga vertical contra el alargamiento se compara con los resultados 
reportados por Aguilar et al. (2011), como se muestra en la Figura 16, donde se 
tienen diferencias hasta del 3.09% y 2.93%, respectivamente. 
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Figura 15. Esquema de referencia: 


7 



b) 

a) puntos de control y b) forma deformada. 


a) b) 



Figura 16. Curva de la carga vertical contra el desplazamiento: a) acortamiento 
entre puntos 1-7 y b) alargamiento entre puntos 10-4. 


Los elementos mecánicos calculados en las dovelas durante la primer etapa de 
carga se muestran en la Figura 17, los cuales no alcanzan los valores umbrales. Por 
otra parte, los elementos mecánicos en las dovelas y en revestimiento secundarlo 
durante la segunda etapa de carga se muestran respectivamente en las Figuras 18 
y 19, los cuales tampoco alcanzan los valores umbrales. 



Figura 17. Elementos mecánicos en las dovelas durante la primera etapa de 
carga: a) momento flexlonante (kN-m), b) fuerza cortante (kN) y c) fuerza normal 

(kN). 
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Figura 18. Elementos mecánicos en las dovelas durante la segunda etapa de 
carga: a) momento flexionante (kN-m), b) fuerza cortante (kN) y c) fuerza normal 

(kN). 



Figura 19. Elementos mecánicos en el revestimiento durante la segunda etapa de 
carga: a) momento flexionante (kN-m), b) fuerza cortante (kN) y c) fuerza normal 

(kN). 


Agullar et al. (2011) realizaron la prueba de carga en condiciones de trabajo, por lo 
que las cargas aplicadas, mostradas en la Figura 14, no producen daño en el 
revestimiento primario ni en el secundarlo. Sin embargo, en este caso de estudio, 
las magnitudes de las cargas aplicadas se incrementaron como se muestra en las 
Figura 20, con la finalidad de estudiar la respuesta estructural de un túnel con 
revestimiento simple y con revestimiento doble en el Intervalo no lineal de material. 
Cada esquema de carga, mostrados en la Figura 20, se aplicaron en el túnel con 
revestimiento simple y con revestimiento doble, respectivamente. 
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Figura 20. Esquema de carga en los modelos numéricos del túnel con: a) 
revestimiento simple y b) doble revestimiento. 


Las dovelas del revestimiento primarlo y del revestimiento de concreto colado en 
sitio se modelaron con el elemento viga-columna gruesa con discontinuidades 
embebidas, la malla de elementos finitos en el programa GID se muestran en la 
Figura 21. Por otra parte, las juntas entre las dovelas se modelaron con elementos 
resorte rotacional, a los que se les asignó el modelo constitutivo propuesto por 
Jansen (1983). Los modelos constitutivos asignados a las dovelas, al revestimiento 
secundarlo y a las juntas longitudinales se proporcionan en el apéndice A2. 




Figura 21. Modelos numéricos utilizados en: a) la primera y b) segunda etapa 
de carga para el túnel con revestimiento doble. 

La carga de ovallzaclón contra la ovallzaclón del anillo con revestimiento simple y 
del anillo con revestimiento doble se muestran en la Figura 22a, donde se observa 
que la capacidad de carga máxima en el anillo con revestimiento doble es 
aproximadamente el doble de magnitud que la capacidad de carga máxima en el 
anillo con doble revestimiento simple. La ovallzaclón en el túnel con revestimiento 
simple es mayor que la del anillo con revestimiento doble, lo que se atribuye a la 
rigidez de los revestimientos. Después del desarrollo de cinco articulaciones 
(articulaciones rotuladas de I a III) en el túnel con revestimiento simple y del 
desarrollo de siete articulaciones (articulaciones rotuladas de 1 a 5) en el túnel con 
revestimiento doble, las curvas de la carga de ovallzaclón contra la ovallzaclón 
ablandan drásticamente, respectivamente. El ablandamiento de estas curvas se 
asocia a la formación progresiva de articulaciones en los revestimientos con el 
Incremento de la carga de ovallzaclón. 
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En el túnel con revestimiento simple, la primera articulación ocurre en la parte 
inferior (I). Posteriormente, dos articulaciones ocurren en la parte superior (II) y las 
últimas dos articulaciones ocurren en los lados (III) como se muestra en la Figura 
22b. Un total de cinco articulaciones de desarrollan en el túnel con revestimiento 
simple como en el revestimiento del túnel Botlek Rallway mostrado en la Figura 11 b. 
Por otra parte, en el túnel con revestimiento doble, la primera articulación ocurre en 
la parte superior del revestimiento primario (1); posteriormente, una articulación 
curre en la parte inferior del revestimiento secundarlo (2), en la parte más gruesa. 
Después, dos articulaciones ocurren en ambos lados del revestimiento primarlo (3). 
A continuación, una articulación ocurre en la parte inferior derecha del revestimiento 
primarlo (5). Más tarde, una articulación ocurre en la parte Inferior del revestimiento 
primarlo (6). Después, dos articulaciones ocurren en ambos lados del revestimiento 
secundarlo (7). Finalmente, una articulación ocurre en la parte superior de 
revestimiento secundarlo como se muestra en la Figura 22c. Un total de once 
articulaciones en el túnel con revestimiento doble, siete articulaciones en el 
revestimiento primarlo y cuatro articulaciones en el revestimiento secundarlo. 



Figura 22. a) Curva de la carga de ovalización contra ovalización y ocurrencia de 
articulaciones en b) túnel con revestimiento simple y c) túnel con revestimiento 

doble. 


5 Conclusiones 

La matriz de rigideces de del elemento finito viga-columna 2D con 
discontinuidades Interiores es simétrico, reduciendo el tiempo de cómputo y 
evitando problemas de inestabilidad numérica, lo que garantiza una formulación 
compacta y computacionalmente eficientemente. Los tres GDL adicionales en el 
nodo interno de este elemento, correspondiente a los saltos de desplazamiento, 
pueden ser condensados estáticamente como secundarlos, reduciendo de nueve a 
seis GDL. Así, la matriz de rigideces se expresa en términos de los grados de 
libertad primarlos en ambos extremos del elemento, manteniendo constante el 
número de GDL externos. 

En el elemento viga -columna con discontinuidades interiores desarrollado, las 
articulaciones pueden ocurrir en cualquier lugar dentro de la longitud del elemento, 
lo cual es una ventaja sobre la forma tradicional de modelar la ocurrencia del daño 
en marcos, colocando resortes con comportamiento no lineal del material en los 
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extremos de cada elemento, donde las articulaciones pueden ocurrir, tal que los 
resultados son dependientes de la malla y los modelos requieren mayor consumo 
de cómputo. El elemento viga-columna gruesa con discontinuidades interiores 
modela adecuadamente discontinuidades de desplazamiento axial, rotación y 
desplazamiento transversal, liberando energía adecuadamente como el desarrollo 
de articulaciones. 

Este elemento fue validado con el modelado de dos túneles con revestimiento 
simple. En el modelo del túnel Botlek Rallway se desarrollaron cinco articulaciones 
dentro de las dovelas, mostrando que la ocurrencia del daño modeladas como 
articulaciones son congruentes con las grietas calculadas con el modelo de 
agrietamiento distribuido realizado por Luttlkholt (2007), pues las grietas y 
articulaciones ocurren en el mismo lugar dentro de los segmentos. Además, se nota 
que, como se esperaba, la pendiente de la curva de carga de ovallzaclón contra la 
ovallzaclón se reduce gradualmente conforme las articulaciones se desarrollan 
dentro de las dovelas. Esta pendiente es casi constante después del desarrollo de 
las últimas dos articulaciones. 

La Interacción del revestimiento primarlo con el revestimiento secundarlo genera un 
Incremento en la carga necesaria para producir la primera articulación en las dovelas 
y una reducción de la ovallzaclón del revestimiento. Se muestra que las curvas de 
capacidad calculadas numéricamente que la contribución del revestimiento 
secundarlo produce una respuesta más rígida del revestimiento y un incremento de 
su resistencia, por lo que, el modelo con doble revestimiento tiene además el doble 
de la carga máxima que el modelo con revestimiento simple. La degradación de 
rigidez está asociada con el desarrollo progresivo de articulaciones en el 
revestimiento con el incremento de la carga de ovallzaclón. 

La interacción del revestimiento primarlo con el revestimiento secundarlo además 
provee una distribución de las fuerzas internas. Cinco articulaciones se 
desarrollaron en el túnel con revestimiento simple y once articulaciones en el túnel 
con doble revestimiento. En el túnel con revestimiento simple, la primera articulación 
ocurre en la parte Inferior y las últimas dos articulaciones ocurrieron en ambos lados 
del revestimiento. Por otra parte, en el túnel con doble revestimiento, la primera 
articulación ocurrió en la parte superior del revestimiento primarlo y la última 
articulación en la parte superior del revestimiento secundarlo. 

En los casos de estudio considerados se calcularon las curvas de la carga de 
ovallzaclón contra la ovallzaclón, así como las magnitudes del inicio de cada 
articulación. El cálculo de las magnitudes de carga es una ventaja sobre otras 
simulaciones numéricas reportadas en la literatura, en la que sólo se proporciona la 
curva de capacidad calculada. 

Es Importante mencionar que en los túneles con revestimiento simple no solo el 
daño fue modelado adecuadamente con la evolución de las articulaciones, sino 
también que la curva calculada de la carga de ovallzaclón contra la ovallzaclón 
fueron congruentes con las curvas experimentales reportadas en la literatura, 
destacando que en la solución obtenida con este elemento no hubo problemas de 
aforamiento de esfuerzos. 

El elemento finito con discontinuidades interiores puede utilizarse no sólo para 
estudiar daño en estructuras de concreto reforzado, pero también en cualquier otro 
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material. Sólo es necesario asignar las propiedades mecánicas correspondientes a 
los modelos constitutivos como la fluencia y fuerzas últimas. El desarrollo de una 
superficie de falla que incluya la interacción del momento flexlonante, fuerza 
cortante y fuerza axial es una mejora al modelo constitutivo para trabajos futuros. 
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A Modelos constitutivos de los elementos estructurales 

A.1 Túnel Botlek Railway 

Los modelos constitutivos de las dovelas del túnel Botlek Railway se calcularon con 
sus propiedades mecánicas de los materiales y con las propiedades mecánicas de 
la sección, incluyendo el acero de refuerzo. La relación constitutiva de las dovelas 
se muestra en la Figura Al y la versión discreta se muestra en la Figura A2. Por otra 
parte, las relaciones constitutivas de las juntas longitudinales se muestran en la 
Figura A3. 
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Figura Al. Modelo constitutivo continúo de las dovelas: a) momento contra 
curvatura, b) fuerza cortante contra deformación por cortante y c) fuerza normal 

contra deformación axial. 



Figura A2. Modelo constitutivo discreto de las dovelas: a) momento contra salto 
de rotación, b) fuerza cortante contra salto de desplazamiento transversal y c) 
fuerza normal contra salto de deformación axial. 
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Figura A3. Modelo constitutivo continúo de la junta longitudinal: a) momento 
contra rotación, b) fuerza cortante contra desplazamiento transversal y c) fuerza 
normal contra deformación axial. 

A.2 Túnel con doble revestimiento 

El modelo constitutivo de las dovelas, del revestimiento secundarlo se muestra 
respectivamente en las Figuras A4 a A6, en el que los efectos del momento 
flexlonante, la fuerza cortante y de la fuerza normas se consideraron desacoplados. 
Estos modelos constitutivos se asignaron al elemento viga columna con 
discontinuidades Interiores. Los modelos constitutivos de la junta longitudinal se 
muestran en la Figura A7, en la cual se asignó a los resortes localizados en los 
extremos de las dovelas y de las dovelas llave. 
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Figura A4. Modelo constitutivo continúo de las dovelas: a) momento contra 
curvatura, b) fuerza cortante contra deformación por cortante y c) fuerza normal 
contra salto de deformación axial. 
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Figura A5. Modelo constitutivo continúo del revestimiento secundario: a) momento 
contra curvatura, b) fuerza cortante contra deformación por cortante y c) fuerza 
normal contra salto de deformación axial. 
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Figura A6. Modelo constitutivo continúo en la parte gruesa del revestimiento 
secundario: a) momento contra curvatura, b) fuerza cortante contra deformación 
por cortante y c) fuerza normal contra salto de deformación axial. 
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Figura A7. Modelo constitutivo continúo de la junta longitudinal: a) momento 
contra rotación, b) fuerza cortante contra desplazamiento transversal y c) fuerza 
normal contra deformación axial. 










